Generic features of the wealth distribution in ideal-gas-hke markets 
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We provide an exact solution to the ideal-gas-like models studied in econophysics to understand 
the microscopic origin of Pareto-law. In these classes of models the key ingredient necessary for 
having a self-organized scale-free steady-state distribution is the trading or collision rule where 
agents or particles save a definite fraction of their wealth or energy and invests the rest for trading. 
Using a Gibbs ensemble approach we could obtain the exact distribution of wealth in this model. 
Moreover we show that in this model (a) good savers are always rich and (b) every agent poor or 
rich invests the same amount for trading. Nonlinear trading rules could alter the generic scenario 
observed here. 
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Wealth and its distribution plays an important role in 
society. Economic policies, natural resources, and human 
psychology are certainly important factors which govern 
the distribution of wealth. However, some features of the 
distribution are independent of the details. As pointed 
out by Pareto a large fraction of the wealth in any 
economy is always owned by a small fraction of the pop- 
ulation and vice versa. His empirical formulation, later 
named as Pareto's law, describes that the distribution of 
wealth w follows a power law P{w) = w~'^ . Recent stud- 
ies 0] of wealth distribution in several countries also con- 
firm that it is indeed the behavior for the " rich" , which 
is only about 3% of the population. The rest follow a 
exponential or Gibbs distribution. An interesting anal- 
ogy Q has been drawn between the economic system 
and a system of ideal gas where particles and their en- 
ergies are modeled as agents and their wealth. Collision 
between particles is similar to trading between agents 
where energy or wealth is neither created nor destroyed; 
it is just redistributed between agents. As pointed out by 
Yakovenko0, such a process obviously generates Gibbs- 
like distributions observed for the majority of the popu- 
lation. 

Then what is the origin of power-law for the rich? 
Chakrabarti and co-workers 0, 5] pointed out that sav- 
ing is an important factor which decides and dictates the 
distribution for the rich. In a generic society agents have 
different opinions and concepts of saving and, accord- 
ingly, each agent saves a fraction of his wealth and invests 
the rest for trading. The available wealth is then shared 
randomly between two interacting agents. These mod- 
els generically predict a power-law distribution of wealth 
with 7 = 2. Later studies indicate that 7 is not truly 
universal and can be changed in certain specific cases. 
A strikingly different distribution of wealth is observed 
in a system of like-minded agents (when saving propen- 
sity is the same for every agent), where it is asymmetric 
and peaked below the average. Numerically it could be 
well fitted to a Gamma distribution however, recent 



studies indicate a discrepancy . The exact form of the 
distribution is still an open question and in this paper we 
refer to it as Gamma-like distribution. 

It is important to make a distinction between wealth 
and money, although in this paper we have used one as 
a synonym of the other. Wealth is usually refered to 
things that have economic utility, e.g., money and prop- 
erty. During an exchange of money for tangible property, 
the wealth of the involved agents does not change. Again, 
the effective value of tangible property in terms of mon- 
etary units changes in time violating the conservation of 
total wealth.^Such realistic features of economy are not 
included in 3] and 0, 0|; however, these simple mini- 
mal models remarkably explain the universal two-class 
feature of the wealth distributions. In particular, nu- 
merical studies by Chatterjee, Chakrabarti, and Manna 
(CCM) P clearly suggest an algebraic distribution of 
wealth for the rich. Later studies (mainly numerical Q 
and a few analytical 3 ) , also reveal that algebraic distri- 
butions are generic for the CCM model and its variants. 
Exact results, however, are far from reach. In this paper 
we will focus on the exact solution of a generic ideal- 
gas-like economy where saving propensities of agents are 
random and distributed arbitrarily. The CCM model is a 
special case where the distribution of saving propensity is 
uniform. Our results clarify why 7 = 2 for most cases and 
also show a way of getting a distributions when 7^2. 
We have also pointed out that in these class of models (a) 
good savers are always rich and (b) every agent invests 
on the average a fixed amount for trading. 

To be more precise about the model let us consider a 
system of L agents having their saving propensity < 
A < 1 distributed as 5(A). The agents are labeled as 
i — 0,1, ... , L — 1 such that Ao < Ai < ■ • • < Xl-i. Let us 
assume that initially total wealth is randomly distributed 
among agents and the average is 



A pair of agents, chosen randomly in this model, ex- 
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change their wealth in the following way. Each agent i 
saves Xi fraction of its wealth and invests (1 — A^) frac- 
tion for exchange. The available wealth for the pair i and 
j is then shared among the agents in a random fashion. 
Thus, 

E, -> A,i?, + (l-6y)[(l-A.)i?, + (l-A,)i?,],(l) 

where < eij < 1 is chosen randomly for each trading 
from a distribution h{e) where the average is r = (cy ). 

It is quite evident that when Xi = 0, the model is iden- 
tical to the ideal gas model, where particles encounter 
random elastic collisions. In this case, irrespective of the 
initial distribution of energy or wealth, in equilibrium 
(after a sufficiently large number of collisions) wealth is 
redistributed according to Gibbs distribution with a tem- 
perature suitably defined hy E. It explains the wealth 
distribution for a majority which follows the exponen- 
tial law. To understand the origin of thepower-law for 
the rich, as explained by the authors of [J| we must in- 
troduce savings. However, note that the particles are 
no more identical once the saving propensity is different 
for different agents. Then, one's saving propensity is his 
identity. It is thus appropriate to study the ensemble of 
such systems. The idea of an ensemble is not new in sta- 
tistical physics. It is a collection of infinitely many mental 
copies of systems prepared under the same macroscopic 
conditions. Physical observables are then needed to be 
averaged over the ensemble to get rid of the microscopic 
fluctuations. Let us study an ensemble of N systems 
labeled by a = 1, . . . , TV, prepared with the same aver- 
age energy E. Thus, in each system {A^} are the same 
whereas the initial wealth is different in different systems. 
Also different is the sequence of pairs {i,j) selected for 
trading and the sharing of available wealth between the 
pair (e,^) during each trading. Thus, it is appropriate 
to find out the distribution of ensemble-averaged wealth 
defined by 
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where E"' is the wealth of an individual i in system a. In 
a large system {L —> go) agent x = i/L has wealth w{x) 
and saving propensity X{x). Since x is uniform in (0, 1), 
X{x) can be obtained from the conservation of probability 
element 



g{X)dX = dx. 



(3) 



Now we can proceed to derive an effective model of 
exchange for wlx). First, note that since the pair of 
agents are chosen randomly for exchange, a specific agent 
X interacts with different agents in different systems in 
the ensemble. Eventually agent x interacts with every 
other agent when N oo. Second, that for a given 



pair X and y, e^y is different in different systems and the 
effective sharing is (exy) = r. Thus, during trading with 
y, the average wealth w{x) of agent x changes to w'{x; y), 
where 

w'ix; y) = [r + X{x){l ~ r)]w{x) + r[l - X{y)]'w{y). (4) 
Since is stationary in the steady state, we must have 



w{x) = / w'{x]y)dy. 



Equation jsj, together with Q), can be solved as 

C 



w{x) 



l~Xix)' 



(5) 
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where C is an arbitrary constant, can be fixed by the 
average wealth E. Explicitly, 



w{x)dx = E. 



(7) 



Note that neither r nor h{e) appears in Eq. 

The distribution of w must satisfy P{w)dw = dx; thus 



P{w) 



dx _ Cg{l - C/w) 
dw 



(8) 



Here the quantity after the last equality is derived using 
Eqs. © and J^Jl. It is clear from ||SJ) that the asymptotic 
wealth distribution for a generic g{X) is P{w) ~ 
with 7 = 2. However one can choose g{X) — g{l — X) 
to get a different power law 7 7^ 2. For example, when 
5(A) — A{1 — A)" defined in the interval < A < 1 one 
gets 



(9) 



which has been reported earlier 

To compare the exact results with the numerical sim- 
ulations, we must understand certain existing numerical 
difficulties of the model. First, that the saving propen- 
sity is never unity. An agent having saving propensity 
A = 1 is a troublesome member of the system who never 
invests and gains by interacting with other members and 
ultimately owns the whole wealth. In usual numerical 
simulations, the maximum saving propensity for a cho- 
sen set {Ai} is g w 1. For example in CCM model 1 — g 
is 0{\/L). We must account for this correction while 
evaluating C. To do so let us look for a generic system 
where the saving propensity is uniformly distributed in 
(0,(7). In this case A(a;) = qx and using the procedure 
described here one gets 



w{x) 



C 



with C 



E 



In(l-g)- 
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FIG. 1: Exact results for the CCM model P) is compared with 
numerical simulations. Practically they are indistinguishable. 
Main figure shows individual wealth of the agents arranged in 
a nondecreasing order of their saving propensity. The inset 
compares cumulative distribution of wealth. Parameters are 
L = lO"* and N = 10*. 



So, for the CCM model C ^Ej ln(L). Alternately one 
can calculate C as follows. The CCM model is equivalent 
to a system of L agents with Ai — i/L. In this case, 
strictly q = 1 - 1/L and thus C = E/\n{L). 

For comparison, we have performed a numerical sim- 
ulation in a ensemble of iV = 10* systems, each having 
L = 10* agents. The saving propensity is chosen from a 
uniform distribution in the interval (0, 1) and is ordered 
in all the systems such that Ao < Ai < . . . Xl-i- Initial 
wealth is also chosen randomly in each system with fixed 
E. Thus the stating wealth of any particular agent is dif- 
ferent in different systems. Average wealth of each agent 
Wi is calculated after 10^ Monte Carlo steps and then the 
distribution P{w) is evaluated. Figure ^ compares w and 
P{w) with exact results. 

The results discussed in Eq. (|10|l are valid for any q, 
where one expects P(w) — C juP' . However, earlier nu- 
merical studies have reported a power-law distribution 
only for large q. The discrepancy can be explained as fol- 
lows. The minimum wealth of the system is i(;(0) = qC 
and the maximum is w(\) — qCjiX ~ q'^). The width of 
the interval q'^C/ (1 — q"^) is quite small for small q unless 
C is large. Since C is proportional to E one must take 
large average wealth to see such a power law. Beacause 
of the finite lower limit in w, the cumulative distribution 
P{'w > m) = 1 — q~^ + C/n shows a power law up to a 
constant. We have plotted P{w > m) + q^^ — 1 vs to in 
logarithmic scale in Fig. |2] which perfectly fits to C/to 
(see figure caption for details). 

It is interesting to ask what happens when two dif- 
ferent economies interact? For example take a sys- 
tem of L agents, of which q and 1 — g fractions be- 
long to two different organizations c = ±. Their sav- 
ing propensity distributed according to go- (A), respec- 
tively, with < A < 1. During trading agent i of type 



FIG. 2: w{x) obtained from simulations of a restricted saving 
model (0 < Ai < q) is plotted together with the exact result 
l|10|l . Inset compares the cumulative distribution of wealth. 
Here g = 0.8, L = 10* = A^, and E = 50. 

a interacts with agent j (type <t') and shares eij{a,a') 
and 1 — £ij{a,a') fraction of "available wealth", respec- 
tively. Let us assume that {eij{a,a')) = raa' (obviously 

r-i = r-+). The distribution of the grand system is 

now g{X) = (7.9+ (A) + pg^{X) and the effective sharing 

is r — q^r++ + (1 — q^r -|- 2g(l — No matter 

what the value of r is, Eq. (0 is still valid and thus the 
distribution of wealth is given by Eq. ©. So the power 
law P{w) ^ w^'^ and even Eq. ^ are quite robust. 

Equation ijHJl is the central result of this paper, which 
states that the wealth of an agent having saving propen- 
sity A is inversely proportional to 1 — X, irrespective of 
what the distribution g(A) and h{e) are. It clearly indi- 
cates that on the average each agent, independent of how 
rich or poor he is, invests a constant wealth C (which is 
of course (1 — A) fraction of his individual wealth) for 
trading jljjj . And then with equal probability he receives 
r or 1 — r fraction of the available wealth 2C. Thus, 
on the average, the individual wealth is preserved in the 
steady state. 

One can instead write Eq. ((HJ as 

wiX) = with = E-' I (11) 

^ ^ 1 - A J I- X ^ ' 

thus, better saving means better wealth. To verify the ro- 
bustness of Eq. (fTT|l let us divide the system of L agents 
with their savings distributed uniformly in (0, q) into two 
groups : (a) the poor savers who have Xi < z and (b) the 
rich savers who have z < X < q. If the poor savers in- 
teract only with the poor and the rich savers interact 
only with the rich, clearly the system breaks up into two 
independent subsystems of size zL and {q — z)L, respec- 
tively. Correspondingly wealth for the poor and the rich 
are Wp{X) = Cp/(1 — A) and ^^(A) = (7^/(1 — A), where 
Cp and Cr are to be determined independently from the 
initial average of wealth in each system. The total distri- 
bution of wealth is P{w) = zPp(w) -f (g — z)Pr{w)^ where 
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FIG. 3: (a) iu(A) for the poor-rich system, where poor savers 
own 2/3 of the total wealth. Inset shows corresponding 
cumulative distribution of wealth, (b) w{X) for the poor- 
intermediate-rich model compared with CCM model Q- Here 



0.4, b = 0.6, q = 0.95 and, L = 10* 



■ N. 



the distribution for the poor and the rich are Pr^p{w) — 
Cr.p/w^, nonzero in the interval Wp(Q) < w < Wp{z), 
and Wr{z) < w < u'p(l), respectively. Depending on the 
choice of parameters one can get different intervals vifhere 
Pp, Pr , both or none of them contribute to P{w). For 
a certain choice, as described in the inset of Fig. |21 it is 
possible to obtain a cumulative wealth distribution which 
resembles the one observed in reality. Correspondingly a 
poor- rich break up at A = 0.4 is shown in Fig. Of a). This 
discontinuity disappears once some intermediate savers 
(whose saving propensities are limited to z < Ai < 6) are 
introduced who can interact with both the poor and the 
rich savers. Surprisingly w(A) and P{w) in this system 
are identical to that of the original CCM model where 
every agent interacts with every other agent. Fig. Otb) 
compares w{x) obtained from a numerical simulation for 
the poor-intermediate-rich system with the exact results. 
This suggests that that each agent in the system stands 
by their own. Irrespective of their distribution of saving 
propensity and interaction with other agents, any tagged 
agent who has saving propensity A acquires wealth which 
is inversely proportional to (1 — A). Equation is thus 
quite robust. 

So far we have discussed the distribution of wealth 
in a ensemble of infinitely many identical systems. In- 
stead, if we look at any given system which is a member 
of the ensemble, wealth of any particular agent (saving 
propensity A) would show fluctuations about the average 
w = C/{1 — A). Such fluctuations have been studied in 
These numerical studies indicate that the distri- 
bution of fluctuations P{E) for any tagged agent (sav- 
ing propensity A) is a Gamma-like distribution which is 
asymmetric about the mean w = C/(l — A) and is peaked 
at Ec < w. It is also known that the distribution becomes 
symmetric with Ec approaching w when A ^ 1. Since 
usually agents are crowded around the peak of the dis- 
tribution, in any particular system in the ensemble P{E) 



is not very different from P{Ec)- Thus for the rich (who 
have A « 1), the distribution of wealth P{E) is the same 
as P{w) (as Ec w w), which explains the Pareto law for 
the rich in every system in a ensemble Whereas 
deviation of P{E) from a power law is expected for the 
poor. An exact study of fluctuations could reveal the 
discrepancy. 

What happens in a system of identical agents (each 
having same saving propensity A)? Note that we need 
not consider identical copies of the systems now. The 
system itself is an ensemble of identical agents. Clearly 
the average wealth is w — C/(l — A) = E, and thus 
the probability distribution P{w) — S{'w — E). Agents 
in the system differ by their fluctuations and thus the 
distribution of wealth at any given instant would only 
count the fluctuations about the average E, which is not 
different from the distribution of fluctuations of a tagged 
agent in CCM model. Extensive numerical simulations 
4] have shown a Gamma-like distribution in this case. 
Further analytic study might shed light on the exact form 
of distribution. 

In conclusion, we have provided an exact solution 
to the ideal-gas-like markets using Gibbs ensemble ap- 
proach. We point out that in a system of nonidentical 
agents, it is appropriate to consider infinitely many iden- 
tical copies of the systems differing by their initial condi- 
tions. Such ensembles represent the evolution of several 
identical systems under the same macroscopic conditions 
and, thus, physical observables make sense only as an 
ensemble-averaged quantity. One single system, instead, 
would encounter fluctuations which are sometimes incal- 
culably complex. A system of agents having the same 
saving propensity is such a case, where average quanti- 
ties like wealth are identical for every agent and fluctua- 
tions are the only things that count. The central result 
revealed from our exact solution is that in ideal-gas-like 
markets, irrespective of the details of the interaction and 
distribution of saving propensity, an agent having saving 
propensity A would, on the average, acquire wealth which 
is inversely proportional to 1 — A, i.e., better savers are 
richer. Thus, every agent, poor ( small A) or rich (large 
A), on the average invests the same amount for trading, 
contrary to the real economy where rich agents usually 
invest more. To alter such scenario one might modify 
this minimal model so that the investment is non-linear 
in w [note that current investment (1 — X)w is linear]. 
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